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Abstract 

We describe an extension of the nonhnear integral equation (NLIE) tehnique 
to = 1 superminimal models perturbed by $13. Along the way, we also 
complete our previous studies of the finite volume spectrum of the N = 1 
super symmetric sine-Gordon model by considering the attractive regime and 
more specifically, breather states. 



1 Introduction 



Finite size effects play an important role in the investigation of quantum field theo- 
ries. They provide a possibility to determine many important physical characteristics 
of the models such as S-matrices, mass ratios, relations between parameters appear- 
ing in the ultraviolet (UV) and infrared (IR) descriptions of the theory and a great 
deal of qualitative information on the spectrum. The investigation of finite size ef- 
fects in integrable 1-1-1 dimensional integrable quantum field theories is particularly 
useful due to the wide range of exact techniques applicable. 

Several exact methods have been worked out to study finite size effects in inte- 
grable quantum field theories [H [21 El S] • Among all these methods the so called 
nonlinear-integral equation (NLIE) technique [31 Hj proved to be the most efficient 
method to study finite size effects in this family of quantum field theories. In most 
of the cases these equations can be derived from an integrable Bethe Ansatz solvable 
lattice regularization of the model. 

With the help of the NLIE technique the best studied model is the sine-Gordon 
model. The integrable lattice regularization starting from which the NLIE can be 
derived is the inhomogeneous 6- vertex model with alternating inhomogeneities [1]. 
The NLIE is derived for the lattice model, then tuning appropriately the lattice 
constant and the inhomogeneity parameter one gets the NLIE for the continuum 
sine-Gordon model [5]. This NLIE was shown to describe the full finite size spec- 
trum of the sine-Gordon model [HI E] , and moreover considering twisted versions of 
the NLIE one can describe finite size effects of massive integrable perturbations of 
Virasoro minimal models (unitary and non-unitary) [8]. 

In this paper we investigate the = 1 supersymmetric sine-Gordon model 
(SSG) by means of the NLIE technique. By analogy with the sine-Gordon model we 
consider the inhomogeneous 19-vertex model [H] with alternating inhomogeneities 
as integrable lattice regularization of the model. Then using the auxiliary functions 
introduced in [THj and the functional relations among them provided by the T-Q 
relations, the NLIE governing the finite size effects of the regularized model can be 
obtained [H]. Carrying out the prescribed continuum limit procedure one arrives 
at the NLIE of the supersymmetric sine-Gordon model (continuum field theory) 
[12]. In [12] the finite size effects in the repulsive regime of the supersymmetric 
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sine-Gordon model has been studied. 

In this paper we complete this earlier analysis by the extension of the NLIE 
technique to the attractive regime of the model as well as to = 1 superminimal 
models perturbed by $13. 

The paper is organized as follows: In section 2 after recalling some basic gener- 
alities about the supersymmetric sine-Gordon model we present the twisted NLIE 
being valid either in the attractive or in the repulsive regime. Section 3 is devoted to 
the infrared analysis of the NLIE reconstructing the breather S-matrices. In section 
4 the UV limit of the finite size spectrum predicted by the NLIE is computed and 
shown in the twistless case to agree with the spectrum of the expected modular 
invariant c = 3/2 conformal field theory (CFT). The relation between the twisted 
NLIE and the = 1 superminimal models perturbed by $13 is discussed in section 
5. Section 6 is devoted to discuss some examples of breather states of the super- 
symmetric sine-Gordon model. In section 7 we perform some numerical checks on 
our twisted NLIEs. Finally our conclusions and perspectives for future work can be 
found in section 8. 



2 NLIE for the supersymmetric sine-Gordon model 
in the attractive regime 

2.1 The model 

The = 1 supersymmetric sine-Gordon model (SSG) is defined by the action: 



AssG= dt dx I -d^'^d^'^ + iij-f'^d^ij + fitpij cos -1^ + —cos f3(p\ , (2.1) 



where (p is a real scalar, ip is a Majoranna fermion field, /3 is the coupling constant 
and the dimensionful parameter fi defines the mass scale, which can be expressed 
by the kink mass Ai \13\ as follows: 





(2.2) 
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For later convenience, we define a new parameter p by 

In the repulsive regime p > 1 while in the attractive < p < 1. The k-th breather 
threshold is p = 1/k. 

The UV limit of the theory is a c = 3/2 CFT being composed of a Neveu-Schwartz 
(NS) and a Ramnond (R) sectors corresponding to antiperiodic and periodic bound- 
ary conditions imposed on the Majoranna fermion respestively. In the NS sector the 
primary fields are given by vertex operators 

V(%\z, -z) = : e4(t+^)^(^)+(t-^)^(^-)] :, (2.4) 

where r, f G {0, 1}, ipo = 1, ipi = ip of conformal dimensions 

A{nr-) = 1(1 + + - A(^'^") = -(l- !I^V + - 

"•"^ 2\r 2 2 ' 2\r 2 ' - 



while in the R sector the vertex operators corresponding to the primary fields have 
conformal weights 

1 / n mR\^ 1 -r 1 /n mR\^ 1 

and read as 

/?„,^(z, -z) = aiz, -z) : e4(t+^)<^W+(f -^)^(^-)] ,^ (3.5) 

with cr(z, z) being the Ising spin field. 

The allowed set of integer or half-integer valued quantum numbers n, m of the 
above vertex operators (12. 4p . (12. 5p are encoded into the modular invariant partition 
function as follows: 

z{R) = ^7^|(xo(g)xi/2(g) + xi/2(g)xo(g)) J] (2.6) 
+ (lxo(g)P + lxi/2(g)r) Yl +\xi/M\' E j^' 

neZ,me2Z 2n-me2Z+lJ 

where 

n,m 2\R 2 



are the conformal weights of the Gaussian part of the CFT, 77(g) is the Dedekind 
function, {xo{(i),Xi/2{q),Xi/16{q)} are the Virasoro characters corresponding to the 
irreducible representations of a c = 1/2 CFT, q = e^'^''^, r being the modular 
parameter. The compactification radius R related to p as follows: 



fi=yi + H, (2.7) 
2.2 Nonlinear integral equations 

We briefly recall the NLIE for the (twisted) SSG. We use the notations and con- 
ventions of the papers [12] and [11], which the reader is invited to consult for more 
details. We put the SSG on a cylindrical spacetime, with infinite time direction and 
compact spatial extension of length (volume) L. The NLIEs governing the finite 
size effects of the (twisted) SSG read as follows: 

log6(^) = a + iD{e)+ig^{e)+igbie) + iG*^\nB){e)-iG*-\nB){e) 

+ \im(K+^-' *\nY)(e) (2.8) 

hgyie) = Cy + igyie) + {K+U^\nB)ie) + iK~U^\nB)ie), 

where Y{9) = l+y{9), B{6) = 1 + 6(6') and B{9) stands for the complex conjugate 
of B{6). We introduced the notation for any function / 

The eqs. (12.81) contain three type of convolutions, one of them is the usual one 
containing integration along the real axis 

00 

if*9)ix)= / dyf{x-y)g{y), 



while the other two ones are defined by integrating on the complex plane along the 
integration contours T{t) and T{t) (t G M): 

(/ *r 9){x) = j f{x - z)g{z), (/ *- g){x) = j_dz f{x - z)g{z), 

where the curve T{t) is the complex conjugate of r(t). The continuous non-self- 
intersecting contour r(t) fulfills the following properties: 
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1. Rer(±oo) = ±00, 



2. < Imr(t) < min(l,p) 7r/2 Wt G M. 



The kernel functions G and K of fl2.8D read as 



We also introduce the odd primitives of the kernel functions (see appendix D. for 
the choice of branch cuts of Xk{^)), 

e e 
X{9) =2n J dx G{x), xk{9) =271 j dx K{x), (2.10) 

that are important in writing the source terms containing information on the exci- 
tations: 
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(2.11) 
(2.12) 



gy{e) = Urn Qy (e + - ir]) , 



Nh Ms 



9yio) = Y^K{e-h,) + J2ixK{e-v,) + XK{e-v,))-J2(xK{e-s,) + XK{e--s,)) 
j=i j=i j=i 

Mc Mw Msc 

- Y.^K{e-c,)-Y,XKn{e-w,)-Y,XKn{e-w^J^), (2.13) 
j=l j=l j=l 

where the second determination of any function: fn{0) is defined as in [U] 

f/(^) + /(^-^7rsign(Ime)) l<p 
fjj(e) = )-'^ ' ■'^ ^ ^ " ^ (2.14) 

/(^) -/(^-^J'7^sign(Ime)) 0<p<l. 



The objects appearing in the source terms of fl2.8l) are as follows 

1. Type I holes: {hf}, j = l,...,Ni 

2. Holes: {hj}, j = 1,...,Nh 

3. " Close source objects" : {cj} j = 1, . . . , Mc, satisfying the condition: 
Imr(ReCj) < |lmcj| < min(l,p)7r. 

4. Wide effective roots: {wj}, j = 1, . . . , Mw 

5. Self- conjugate effective roots: {w^J}}, j = 1, . . . , Mgc 

6. Ordinary special objects: {sj}, and their complex conjugates {sj} j = 1, . . . , Ns 

7. Virtual special objects: {vj}, and their complex conjugates {vj} j = 1, . . . , Ny. 

Using the light-cone lattice approach [4J the NLIE fl2.8p was derived from the in- 
homogeneous 19-vertex model with alternating inhomogeneities. Thus the source 
objects of (12.81) are related to the Bethe roots and the zeroes of transfer matrices of 
the fusion hierarchy of the underlying 19-vertex model. This relation and the de- 
tailed description of ordinary and virtual special objects can be found in ^LA\. The 
driving term bulk contribution in the equation for In b{6) reads as 

D{e) = i sinh^, 

where the dimensionless scale parameter £ = AiL with Ai and L being the kink 
mass and the volume respectively. 

The values of the constants of the NLIE (12.81) are as follows 

Cb = iinSb + a), a = uj(l + -] +Xoo{N^-N+), Sb^ {0,1}, 



p 

I 3S Scji r s 



p + 2 TT 



S UJ 

T - 



, TT / 1 

Xoo = x{+oo) = - 1 - - 
2 V P 



_P + 2 TT 

and 

Cy = in{S + 6y) + inQ{p- 1) {Mw + M^c), Sy e {0, 1}, 
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where [...] stands for integer part, Q{x) denotes the Heaviside function and u is 
the twist angle of the underlying solvable lattice model. For conventions used here 
see section 3 of ref. [12]. In addition to fl2.8p we need two other equations for the 
determination of type I holes and effective wide and self-conjugate roots. 

The positions of type I holes can be determined from the function a{6) given by: 

-loga(^) = i6a{e) + {K*^\nB){e)-{K*-\nB){e)-Cy, < |Im^| < ^, (2.15) 

Nh Ms 

Sa{e) = J2^K{e-h,) + J2ixK{e-v,) + XK{e-v,))-Y,ixK{e-s,) + XK{e--s,)) 
j=i j=i j=i 

Mc Mw Msc 

- Y.^K{e-c,)-Y,XKn{e-w,)-Y,XKiAe-w'f^). (2.16) 
j=i j=i j=i 

The function necessary to know for the determination of wide and self-conjugate 
effective roots is as follows, for min(l,p) tt < Im^ < ^I^^^tli: 

loga(e) = tDjj{e) + igui{e) + igMi{e) + {Gjj*,lnB){e) - {Gii *^\nB){e) 
+ ((ir-t),,*lnr)(^^) + a, (2.17) 

where 

C^^h'-^'-(^'~"*) ,2.18) 

[iTT{N_-N+) 0<p<l. 

The source objects appearing in the NLIE are not arbitrary parameters, but they 
have to satisfy certain quantization conditions [TT|, [T21 [H] . These are as follows: 

• For holes: 

^ \ogbih,) = 27r4^,, J = 1, ...,Nh. (2.19) 

• For ordinary special objects: 

Imlog6(s,) = 27r/,^, \bis,)\ > I, (Imlog6)'(s,) < 0, j = l,...,Ns. 

(2.20) 

• For virtual special objects: 

lm\ogb{v,) = 2nh^, \biv,)\ > I, (Imlog6)'(t;,) > 1, j = l,...,N^. 

(2.21) 



7 



• For close source objects (only for the upper part of the close effective pair): 

- log6(c]) = 27r/^T, 3 = l,...,Mc/2. (2.22) 

• For wide effective roots: 

- \ogd{w]) = 27r/^T, J = 1, ...,Mw/2. (2.23) 

• For self-conjugate effective roots: 

- \og~a{w]i^^) = 27rJ^T(.), 3 = 1,...,M,,. (2.24) 

So far we have determined only the upper part of the complex pairs, but the other 
parts can be determined by simple complex conjugation. 



Finally for type I holes: 

- loga(/if)) =27r/^a), J = l,...,iVi. (2.25) 



All the above quantum numbers Jq^s are half integers. A state is then identified by 
a choice of the quantum numbers {Ih^ Icp ■■■)■ We also mention that the NLIE itself 
can impose constraints on the allowed values of some of these quantum numbers. 
The counting equations in the continuum theory read as: 

Nh + 2N^ - 2Ns = 2S + Mc + 2e{p - l){Mw + M.J, (2.26) 

Ni - 2N^ = S + Mi-MR-6y, (2.27) 

where S is the topological charge of the theor}l| , Mi denotes the number of effective 
roots (or source objects) with the property | < |Im6'| < ZLi^t^^ jvf/j stands for the 
number of those real zeroes of the function 1 + a{6) which are not type I holes. In 
the lattice Bethe Ansatz they correspond to real Bethe roots. Finally denotes 
the number of real special objects. They are such real objects where a{6) of f l2.15p 
takes the value —1, and i ^ In 0(6*) is negative. 



"'^In this convention charged kinks carry topological charge ±i. In our case S can take only 
integer values because in finite volume the periodic periodic boundary conditions restrict the 
number of kinks to be even. 
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In the repulsive regime, the complex roots describe the internal degrees of free- 
dom of kinks. In the attractive regime, however, configurations consisting entirely 
of wide and self-conjugate effective roots describe the breathers and it is clear from 
fl2.26p that they do not contribute to the topological charge. 

The energy and momentum of a state can be expressed by the solution of the 
NLIE as: 

/ Nh nS 

E = ^A \ cosh(/ij) + {cosh(t>j) + cosh.{vj)} — {cosh(sj) + cosh(sj)} 

Mc Mw Mw 

- Y cosh(9) - Y cosh{w j)ii - Y cosh{w'^J)ii (2.28) 
j=i j=i j=i 



— [ desmh{9) \nB{e) - — / desmh{e) lnB{e) 
2tt Jr 27r Jf 



(Nh n^ Ms 

sinh(/ij) H — ^ {sinh(sj) + sinh(sj)} — ^ {sinh(sj) + sinh(sj)} 
j=i j=i j=i 

Mc Mw Mw 

- Y sinh(9) - Y smh{wj)n - ^ sinh(w^{))// (2.29) 
j=i j=i j=i 

+ ^ j d9 cosh{9) \nB{9) - ^ j_ d9 cosh{9) In B{9) 
We notice that there is no bulk energy term due to the supersymmetry of the model. 



3 Infrared limit and breather S-matrices 

In the infrared limit (L —>■ oo) the bulk driving term D{9) develops a large imaginary 
part giving restrictions on the imaginary parts of the positions of close source objects 
in the infrared limit. In the repulsive regime in (12.171) Djj{9) = thus there are no 
restrictions on the positions of the wide effective roots, only the close source objects 
are forced to form quartets or 2-strings p^. However in the attractive regime Dii{9) 
in 02.171) is no longer zero imposing constraints on the imaginary parts of the wide 
effective roots as well. Analyzing the NLIE and following the argumentation of [6], 
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the effective roots fall into the configurations as follows: 



1. Arrays of the first kind are effective root configurations containing close source 
objects and wide effective roots as well. These type of arrays consist of effective 
roots of the form: 

6k^^ = 6±i{^ — kpn) , 9^^ = 9 ± i {n — fi — {k — 1)ptt) , k = 0, . . 



2p 

with 9 and < /i being real parameters. At certain special values of /i these arrays 
degenerate. There are two degenerate cases: odd degenerate arrays, which contain 
a self-conjugate effective root at 

and accompanying complex pairs at 

.Tt{l - {2k + l)p) 



9. = 9±i- 



k = 0,..., 



2p 



and even degenerate ones, which contain complex-pairs of effective roots at the 
positions 

~ „ 7r(l — 2kp) 
9k = 9±i^ ^, A; = 0, ... 



2P. 



These degenerate arrays contain exactly one pair of close source objects. According 
to the counting equation f l2.26p the presence of close source objects entails the ap- 
pearance of holes which correspond to the rapidities of the kinks. Thus these arrays 
of the first kind describe multiparticle configurations containing kinks and antikinks. 



2. Arrays of the second kind contain only wide- and self-conjugate effective roots. 
They describe the breather degrees of freedom of the model. There are two types of 
them. The odd ones contain a self-conjugate effective root 

and wide effective-pairs at 

^~. = ^±.^^il^i^±M, A; = 0,...,., 0<.< 
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while the even ones contain only wide effective-pairs 



7r(l - 2kp) 
2 



A; = 0, 0<s< 



2p 



1 



1. 



They correspond to the (2s + l)-th breather -B2S+1 and the (2s + 2)-th breather 
-B2S+2, respectively. 

The deviations of the imaginary parts of the previous effective root configurations 
from the formulae listed above are exponentially small in the volume L. 

It can be seen that arrays of the second kind become degenerate ones of the first 
kind, if one analytically continues increasing p. The reason is that breathers are of 
course kink-antikink bound states, while degenerate arrays of the first kind describe 
scattering states of a kink and an antikink. 

In the infrared limit one can drop all terms containing the convolution of \aB{9) 
and In B{6) , because they are exponentially small in the volume L. One can therefore 
compute the energy and momentum contribution of an array of the second kind 
corresponding to the breather Bg. The energy- momentum contribution turns out to 
be 



where 9 is the common real part of the roots composing the array. This is just the 
contribution of a breather B^ moving with rapidity 9. Arrays of the first kind do 
not contribute to the energy and momentum in the infrared limit lending support 
to their interpretation as polarization states of the kinks. 

Now we proceed to show that with the above interpretation the NLIE correctly 
reproduces the 2-body scattering matrices of the SSG including breathers. In the 
repulsive regime the 2-kink S-matrices have been calculated in detail in |12j . 

Hereafter we will concentrate only on breather-breather 2-body scatterings, be- 
cause from the NLIE or equivalently from the study of the finite size effects of our 
model the direct determination of kink-breather 2-particle scattering amplitudes is 
not possible. This is because according to the counting equation (12.261) only even 
number of kinks and antikinks are allowed in finite volume. So the simplest scat- 
tering amplitude containing breathers and kinks too, which can be extracted from 
the finite volume analysis of the model contains 2-kinks and one breather, which is 
a 3-body scattering amplitude. 
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3.1 Breather S-matrices 

For the sake of simplicity we will show that the NLIE correctly reproduces all the 4 
eigenvalues of the 2-body scattering matrix of the first breathers. In the language of 
the NLIE the first breathers are represented simply by self-conjugate effective roots. 
Let the two self-conjugate roots corresponding to the two Bi breathers: 

^« = ^i + ^^^^^, and w^^J = e, + z^^^P±^, e,,e,eR. 

Then the S-matrix eigenvalues can be read off from the quantization conditions 
of 6i and 62 given by the auxiliary function loga(^). In the infrared limit, as we 
have already mentioned before, all terms containing the convolution of lni?(^^) or 
In B{6) can be dropped, nevertheless the convolutions containing Iny(^) must be 
kept, because Y{6) tends to a finite non-trivial function in the infrared limit. So in 
the infrared limit the quantization condition (12.241) reads as 

loga(«;(^)) = tismhniw^^J) + 4°Vic^) - ^iXn^M^ - ^i? - 

+ z J2 XKiiiw^^ - h^) - ^{x^I%I{-^l), (3.1) 

where 

00 

= jdy {K--/')n{wi'J - y) \nY{e). (3.2) 

—00 

The bulk source term of (13. ip gives the phase shift of the two breather wave function: 

i£sinh//(w£-') — > —i 2Ai sin — L sinh 6*1, 

while the third term on the right hand side of (13.11) always provides the sine-Gordon 
part of the expected S-matrix eigenvalue: 

(See appendix A. for the SSG breather-breather S-matrix amplitudes.) The other 
terms in (13. ip describe the supersymmetry factors of the breather S-matrices. The 
four eigenvalues of the Bi — Bi S-matrix correspond to states with different values of 
5y and different number of type I holes. According to the counting equation (I2.27P 
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in case of 6y = (Neveu-Schwartz sector) the value of A^^i can be either zero or 2, 
while in case of 6y = 1 (Ramond sector) A^^i is either zero or 1. 



The 6y = 1 case: 



In this case in the infinite volume limit the function Y{9) reads as 
Y(^y=^)(e) = sin (f ) cosh {'-^) cosh(^ - 9,,) ^9,-9, 



n^sinh + If - ^) sinh + ^ + ^) 



2 

(3.3) 



Then substituting this into the integral if"* {wgc) (13.21) one gets: 



(0) _ cosh (g^) , ^ cosh (g^) , cosh (g^gipPE) 

h [Wsc)\^^^^^,_,,^ - ^^^1^ + in ^^^j^ cosh (g^ME^) 

- ^XG{0-0l)-^XG{O-92), (3.4) 

where 



oo 



Xc(^) - / - sm(^g) —^^^^^^y^^^^^^^- (3-5) 







Then inserting this result into (13.11) and taking its exponential in the A^i = case 
one gets the eigenvalue A^^(^^i — 92) of flA.Sp of appendix A. In the A^^i = 1 case the 



position of the single type I hole in the infrared limit is h^i^ = 9i2. Then putting ev- 
erything together the quantization condition (13. ip gives the eigenvalue A^_}\9i — 92) 
of (ESD. 



The Sy = case: 



In this case in the infinite volume limit the function Y{9) reads as 

(,^^0^ cosh (2(9 - 9,2)) + cosh(gi -92) + 2 sin' (f ) ^1-^2 



4 n^sinh + ^ - ^) sinh + ^ + ^) 



(3.6) 

Then substituting this into the integral I^\wsc) (13. 2p one gets the result summa- 
rized in appendix B. Then inserting the result of this integral into the quantization 
equation (13.10 and putting everything together one gets that in case of A^i = the 
S-matrix eigenvalue A^^\9i — ^2) of (1A.9P is reproduced correctly, while the other 
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eigenvalue A^''(^i — ^2) of ( ]A.9p is obtained from the Ni = 2 case where the positions 
of the type I holes are h^^ = 612 + \nA± with A± is given by (lB.7p . 

Thus we demonstrated that all the 4 eigenvalues of the Bi — Bi S-matrix can 
be reproduced by the NLIEs f l2.8l) -( l2l25l) . All the 4 two Bi states are described by 
the same effective root configuration, but they differ in the number of type I holes 
and in the choice for 6y. To close this section we also mention that the breather 
S-matrices can be determined solely from the function loga(^), so the breather S- 
matrices calculated from the NLIEs are insensible for the value of the constant 6b 
and a of (12.81) . However, as we will see in the next section, these parameters strongly 
influence the UV conformal weights of the states. 



4 UV limit and kink approximation 

In this section in the presence of a twist angle u we examine the UV limit of the states 
described by the NLIEs ( I2.8l) - fl2.25l) . We follow the standard approach described in 
detail in [3l [6l [261 [2Z]- The positions of the sources for £ — > can remain finite 
(central objects), or they can move towards the two infinities as ±ln (|) (left/right 
movers). We introduce the finite parts of their positions 6j by subtracting the 
divergent contribution: 

W}-{9f±ta(^).9»}. 

We denote the number of right/left moving and central objects by N^'^, Np^, M^'^, . , 
etc. For later convenience we introduce the right/left moving and central spin given 
by 

^±'0 = l(Ar±'0 + 2ArJ^'° - 2Ar^'° - M^'° - 2 0(p - 1)(M^'° + M^f)). (4.1) 

In the UV limit the NLIE splits into three separate equations corresponding to the 
three asymptotic regions. This is why for all the auxiliary functions of the NLIE 
(12.81) we define the so called kink functions as 

F±{e) = \imF (e ±\n'^] . F G {logfo, logy, logo, . . . }. (4.2) 
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In the UV limit these kink functions satisfy the so called kink equations which can 
be straightforwardly determined from f l2.8p . The energy and the momentum in the 
conformal limit can be expressed by the left and right kink functions. 

Following the method worked out in [3l El |27] it turns out that in the UV limit 
the energy and momentum can be expressed by a sum of dilogarithm functions with 
the 6 ±00 limiting values of the kink functions in their argument. One group of 
these two limiting values is trivial, namely: 

6±(±oo) = 6±(±oo) = 0, (4.3) 

l/±(±oo) = (-1)^^ (4.4) 

The other limiting values satisfy a nontrivial coupled nonlinear algebraic equations, 
the so-called plateau equations [B]: 

±^ _L -U ( l\/r I ^1 AT± 



log 6± (too) = a ± 2^X00 (5 - 25±) ± 2m ± m (^M,, + ^ - Nl 

+ ^ {ln5±(Too)-ln5±(Too)} + ilny±(Too), (4.5) 

2/±(Too) = [-Ifs^+'y (5±(too)5±(too))^ (4.6) 

where are integers depending on the relative positions of the complex effective 
roots and their actual value is irrelevant from the point of view of the exponent of 
equation (14. 5p . 

To solve eqs. fl4.5|4.6p analytically we take the Ansatz of ref. 



h±{Too) = e^^'P^ 2 cos(p±), 6±(too) = e^^^^± 2 cos(p±), (4.7) 

B^iToo) = e-- B^iToo) = e^-- (4.8) 

sm (p±) sm (p±) 

y±{Too) = r±(Too) = 4cos(p±)^ > 0. (4.9) 

sm (p±) 

Since we need only b±{^oo) and not its logarithm, we can restrict the allowed values 
of p± in the [0, 27r] interval. The solutions of the plateau equations formally take 
the form 

p± = n(kf±5b±- + A,^ - (2kf ± 25b ±2- + 3Ap^ , (4.10) 
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where 



A,^=S-2S^~Np, 







rp±i 


= 

) "'p 






rp±i 






. vr . 




. TC . 




. vr - 



• (4.11) 



and further constraints must be satisfied by the parity of the integers and N^: 

(4.12) 
(4.13) 



= 2/J + M,, + ^ - Nt - n% If e Z, 



mod 2 = 2S^ + 6y mod 2 

Using the dilogarithm sum rule of appendix C and putting everything together the 
conformal weights take the form 

2 



I fn±-^ S 
2" 



where 



R 



6h±k^± -A 



N. 



± 



^\ T S^n± - ^{S - S^)S^, = mod 2. 



(4.14) 

(4.15) 
(4.16) 



J± = ±4± ± 24± =F 2/s± =F 2/c±T =F 2/^±T =F / ±t ± IhW± 



T (S^ + e(p - 1)(M± + M,t)) ('^f> ± ^ T ± M,, T (4.17) 
T (^)%^^.T(M^ + M,t) (- 



where 



M, 



4± = ^ 4± , -?'s± = ^ -''si , -^fi = -^t;± ' 4±T = /^±T, ■ . .etc., 



and K± are integers depending on the concrete complex effective root configuration 
under consideration. Analysing the formulae above it can be proven that there is a 
relation between n± of fl4.15p and the topological charge 5, namely in the 5y = Q 
case: 

n± e Z if ^ G 2Z, (4.18) 
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if 



5 G 2Z+ 1, 



(4.19) 



while in the 5, 



y ~ 



1 case: 



n± G Z + - 



if 



S e2Z 



(4.20) 



n± G Z 



if 



5 G 2Z+ 1. 



(4.21) 



Moreover it can be proven that depending on the state under consideration the sum 
N± + J± can be either integer or half-integer, but in the Sy = 1 case the sum N± + J± 
is always an integer. 

Thus, just like in the analysis of the UV conformal weights (14. 14114. iTj) in the 
twistless case {a = u = 0) yields that the operator content of the UV limit of the 
theory with ordinary periodic boundary conditions is that of the modular invariant 
c = 3/2 CFT with partition function given by (12. 6p . The topological charge S of the 
model can be identified with the winding number m of (12. 6p . and 6y is the parameter 
of the NLIE which distinguishes the Neveu-Schwartz (NS) and Ramond (R) sectors 
of the theory, namely Sy = corresponds to the Neveu-Schwartz sector, and Sy = 1 
describes the Ramond sector of the model. 

5 Twisted NLIE and = 1 superminimal models 
perturbed by <l>i3 

Based on earlier results [HI [ISl [HI [IHl E] it is expected that the NLIE (gSD at 
appropriate values of the twist angle a can describe the finite volume spectrum of 
the supersymmetry preserving perturbation $13 of the = 1 superminimal models. 
In this section we confirm this expectation. 

The N = 1 superminimal models SAi{r, s) can be characterized by two positive 
integers r and s such that r and are coprime integers. Their central charge is 



The highest weights of primary fields are also characterized by two integers, k and 



given by 




(5.1) 
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I: 

8rs 32 

(5.2) 

The sum A; + / is even in the Neveu-Schwartz sector and odd in the Ramond sector. 
The special choice s = r + 2 corresponds to the series of unitary superminimal 
models. 

The parameter a of the NLIE is related to the twist uj of the underlying Bethe 
Ansatz solvable lattice model as: 

a = uj^ + Xoo{N^-N+)- (5.3) 
P 

The twisted lattice Bethe equations (see section 3 of ref. [121) are invariant under 
the change u ^ u + n. This discrete symmetry of the integrable lattice model is 
reflected by our NLIE in such a way that when we replace u hj uj + n, according 
to (15.31) the value of a changes by 27r times an integer, but shifting a by 27r is an 
invariance of the NLIE (12.81) . with appropriate redefinition of the auxiliary functions 
(log 6, logy, ...etc.) and the Bethe quantum numbers. 

From now on we restrict ourselves to the case of neutral (i.e. S = 0) states. Then 
the superminimal models SAi{r, s) perturbed by $13 can be obtained by an RSOS 
restriction of the SSG model [21j. This RSOS restriction is similar to the one in the 
ordinary sine-Gordon case [23] • To get access to the perturbed superminimal 
models we should restrict ourselves to the S = neutral charge sector of the theory 
in such a way that the parameter p of the NLIE must be 

2r 

p = (5.4) 

s — r 

and the twist must take the values: 

vr 

UJ = k mod vr. A; G Z. (5.5) 

p + 2 ^ ' 

All of these values of the twist uj is necessary to get all the states of the corresponding 
RSOS model. On the other hand not all these twist values correspond to inequivalent 
values of a and so to different physical states. This is a consequence of the RSOS 
truncation. 
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Following from fl5.3p . fl5.4l) and (15.51) the values of a corresponding to the neces- 
sary RSOS restriction are as follows: 

ak = k—, fc = 1,2, . . . ,r - 1, (5.6) 
r 

where we introduced the notation 



t 



s — r 



2 

The t = 1 special choice corresponds to the case of the unitary superminimal models. 

Restricting ourselves to the S = neutral sector, inserting (15. 6p . (15. 4p and (12. 7p 
into the general formulae of UV conformal weights (14. 14114. l7|) . and correcting the 
final formula by a trivial contribution coming from the difference of the central 
charges of the c = 3/2 and the superminimal CFTs, one gets the UV conformal 
weights as follows: 

A-(„„ i) = ^ + ((2n. + t)r-fa)^-(r-.)^ + iV, + J,. (5.7) 

lb ors 

Exploiting (I4.18P and (14.201) one can see that this formula agrees with (15. 2 p if one 
identifies the integer / of (15. 2p with 2n± + k in i \5.7\) . So as not to overflow the Kac 
table of the superminimal models one should impose the constraint: 

1 < 2n± + A; < s - 1, A; = 1, 2, . . . , r - 1. (5.8) 

Furthermore locality requirements (i.e. A+ — must be integer) force us to ac- 
cept only those solutions of the plateau equations where n+ = Restricting our 
attention to states satisfying the previous two constraints we get the UV conformal 
weights of the SAi{r, s) superminimal models. This fact suggests us that the twisted 
NLIE at the values of the parameter a given by (15. 6p describes the finite size effects 
of the integrable perturbations of the SA4{r, s) models. 

To close this section we remark that the parameter a does not occur in the equa- 
tion of loga(^) in the attractive regime. Since loga(^) is the only counting function 
which determines the breather-breather S-matrices, the infrared asymptotics of the 
breather states does not depend on a and so the S-matrices involving only breathers 
are also unchanged. As a consequence in accordance with [2jL|, the NLIE yields that 
in the perturbations of non-unitary superminimal models, the breather-breather S- 
matrices can be calculated simply by inserting the appropriate value of p from (15. 4p 
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into the formula of the breather S-matrices of appendix A. For example in the second 
section we demonstrated that the Bi — Bi S-matrix can be reproduced by the NLIE. 
According to fl5.4l) the p = 2/3 choice corresponds to the supersymmetric Lee- Yang 
model, so substituting p = 2/3 into the Bi — Bi S-matrix, we can reproduce the 
well known S-matrix of the supersymmetric Lee- Yang model. 

By this simple example we demonstrated that not only the UV limit of the 
perturbed superminimal models are reproduced correctly, but their S-matrix as well. 
Here we remark that in the unitary perturbed superminimal models, the actual value 
of the twist-like parameter a plays an important role in the determination of the 
S-matrix. In this case we have only kink S-matrices and after lengthy calculations 
it can be shown that the twisted NLIE reproduces the kink-kink S-matrices [21] of 
the perturbed unitary superminimal models. 



6 Some examples of breather states 

In this section we consider some of the simplest Bi breather states and with the help 
of the NLIE we determine their UV conformal weights. In the attractive regime 
the simplest excitations of the supersymmetric sine-Gordon model are the single 
Bi breather states at rest. Due to the zero momentum condition and the left- 
right symmetry of the counting functions, the single self-conjugate effective root 
characterizing the excitation is located exactly on the imaginary axis, i.e. Wgc = 
i^{p + 1). Thus no quantization condition is necessary to be imposed on the self- 
conjugate effective root. 

First let us consider the static Bi breather state quantized with 6y = 6y = 0. 
Calculating the UV conformal dimensions along the lines of section 4 we obtain the 
conformal dimensions: 




so the UV limit of this state is the linear combination |1— ) = -^{Vi q (0,0) — 
V"ii'o^(0, 0)|0) of the vertex operators V^i^J of the modular invariant c = 3/2 CFT 

(ED. 

Using conformal perturbation theory (PCFT) the leading order correction to this 
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fe|i_)m (NLIE) 


fe|i_)m (PCFT) 




0.1 


0.028726870911825 


0.02873952098276545 


6.3 • 10"^ 


0.05 


0.009479129734511 


0.00948050632214565 


6.86 ■ 10-5 


0.01 


0.000721896146209 


0.00072190412804693 


3.98-10"^ 


0.005 


0.000238138683370 


0.00023813955194438 


4.33 • 10"^ 


0.001 


0.000018133406805 


0.00001813341184091 


2.51 ■ 10-1° 


0.0005 


0.000005981794545 


0.00000598179509334 


2.73 ■ 10-11 


0.0001 


0.00000045549071 


0.00000045549071160178 


1.58 ■ 10-13 



Table 1: Numerical comparison of PCFT with NLIE for the state |1-) at i?^ = 5. 
state has been calculated in 



and it reads as: 



^E|i„)(£)+(^-12(A+ + A-) 



-d^-^ e^y + o{i^y), (6.1) 



where y 



where 



and the leading order coefficient takes the form: 



a = -. 8^/^' - ( ^ 

2 y(i_^)U^-l 



7(x) 



r(x) 



rri -x) 



and 



^|i-> =7(2 



2i?2 



^ ' i?2 P V 2 + 2i?2 



2i?2 



7 



1 



(6.2) 



(6.3) 



(6.4) 

We checked numerically our NLIE against PCFT at a lot of values of the com- 
pactification radius R, and in every case we experienced very good agreement. To 
illustrate the agreement between NLIE and PCFT, the numerical comparison of 
fe|i_)(£) coming from NLIE and leading order PCFT at the specific = 5 point 
can be found in table 1. One can see that the two sets of data approach to one 
another as the volume tends to zero. 

Another important example of the one Bi breather states is the static one Bi 
breather state in the Ramond sector. This state is quantized by 6y = 1, 6b = 0. The 
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UV conformal weight corresponding to this state is 

= — + Ar, (6.5) 

which turns out to be the hnear combination of the vertex operators R^i o(0, 0). 

To close this section we also calculate the UV conformal weights of some of the 
two Bi breather states with lowest energy and zero momentum. These states contain 
two self-conjugate effective roots. One of them is a left mover and the other is a right 
mover in the UV limit. According to section 3 these states may contain type I holes 
as well. So in this case not only the quantum numbers 6b and Sy characterize the 
lowest lying state but the number of type I holes as well. Performing the calculation 
of the UV conformal weights of the lowest lying zero momentum two Bi states one 
gets that the state having quantum numbers 6y = 6b = Ni = corresponds to the 
symmetric first level descendent of the vacuum with weights: 

A± = l. 

The state quantized with quantum numbers 6y = 6b = and Ni = 2 has UV weights 

A± = 1/2, 

corresponding to the state '?/'-?/'(0, 0)|0). 

In the Ramond sector (i.e. 6y = 1) the two Bi state, which has quantum numbers 
6b = 1 and A^^i = in the infrared limit, will have quantum numbers A^^i = 2, = 1 
in the UV limit due to the fact that the counting function i loga(^) will have a 
negative slope part creating new type I holes. The calculation of the UV conformal 
weights of this state yields: 

~ 16 ^ 8^' 

which turns out to be the same as that of the one Bi state at rest of the Ramond 
sector (16.51) . so this state must originate from the other linearly independent com- 
bination of the vertex operators -R±i/2,o(0, 0) in the UV. 

7 Numerical test of the twisted NLIE 

The most obvious test of the twisted NLIEs for the perturbed superminimal models 
would be the numerical comparison against data coming from truncated conformal 
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space approach (TCSA). Unfortunately there is no available TCSA data in the liter- 
ature for this class of models. This is why we choose another way to test the twisted 
NLIEs. It is well known that the c = 7/10 unitary CFT is present in the unitary 
series of superminimal and Virasoro minimal models as well. In the superminimal 
language it corresponds to SAi{3,5) whilst in the minimal CFT language it cor- 
responds to A^(4, 5) and certainly their integrability preserving perturbations are 
also equivalent. In ref. [8] the NLIE description of the finite size effects in minimal 
models perturbed by $13 was given introducing an appropriate twist parameter into 
the NLIE of the sine-Gordon model. For later convenience we will refer to these 
equations as twisted Destri-de Vega (DDV) equations. So the finite size spectrum 
of the A^(4, 5) + $13 model can be described by two different sets of NLIEs. First 
considering the model as a perturbed superminimal model the set of NLIEs is given 
by the equations (12.8112. 251) with p = 3 and twist values = Then formulating 
the model as a perturbed minimal model the finite size effects are governed by the 
twisted DDV equations of ref. ^8i| with Pddv = 4 and twist values a^^ljy = ^■ 

Hereinafter we will compare numerically the two types of NLIEs for the ground 
states of the SM{3, 5) + $13 model. It is well known that the SM{3, 5) + $13 model 
has 3 ground states which are degenerate in infinite volume, but they split as the 
volume is decreased. In the UV the 3 ground states of the model correspond to 
primaries with conformal weights ^ {0' il)}- language of the twisted 

DDV equations they correspond to twist values: ctDDyG {f,f,x} respectively. In 
the language of our twisted NLIEs (12. 8112. 25]) the ground states with UV conformal 
weights and ^ can be found in the NS sector of the model {6y = 0) and they 
correspond to twist values: a = f and a = ^ respectively. The ground state with 
A^*^ = ^ is in the R sector {6y = 1) and the corresponding twist value is a = 
The numerical comparison of the two types of NLIEs for the 3 ground states of the 
model can be found in tables 2,3 and 4. The high accuracy agreement between the 
two sets of numerical data reassures the correctness of our twisted NLIEs for the 
SM{r, s) + $13 models. 



23 





E^ fNLIE) 


E^ fDDV) 


2 


-0 0597271103 


-0 0597271103 


1 


-0.2399086782 


-0.2399086782 


0.5 


-0.6282004195 


-0.6282004195 


0.1 


-3.6166077208 


-3.6166077208 


0.05 


-7.2973738367 


-7.2973738367 


0.01 


-36.6389894487 


-36.6389894489 



Table 2: The first ground state of the superminimal model <SA^(3,5) perturbed by $13 corre- 
sponding to A='= = 0. 



e 


E2 (NLIE) 


E2 (DDV) 


2 


0.0007298114 


0.0007298114 


1 


0.0109172291 


0.0109172291 


0.5 


0.0613220318 


0.0613220318 


0.1 


0.7610380973 


0.7610380973 


0.05 


1.7468225711 


1.7468225711 


0.01 


9.9655028929 


9.9655028911 



Table 3: The second ground state of the superminimal model iSA^(3, 5) perturbed by $13 corre- 
sponding to = 



e 


E3 (NLIE) 


E3 (DDV) 


2.5 


0.0344866621 


0.0344866621 


2 


0.0668483893 


0.0668483893 


1.5 


0.1380515377 


0.1380515377 


1 


0.3217601530 


0.3217601530 


0.8 


0.4812047615 


0.4812047615 


0.5 


1.0154869429 


1.0154869429 



Table 4: The third ground state of the superminimal model SA4{3,5) perturbed by <I>i3 corre- 
sponding to A='= = jq. For i < 0.5 special objects appear in the NLIEs obstacling the convergence 
of the usual numerical iteration. 
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8 Summary and perspectives 



In this paper we presented and analysed the NLIE governing the finite size effects 
in both the attractive and repulsive regimes of the SSG. We also demonstrated that 
the twisted version of the NLIE at certain particular values of the twist parameter 
describes the finite size effects in the = 1 superminimal models perturbed by $13. 

As a starting point, we extended our previous studies [T2] of the SSG to the 
attractive regime in order to describe breathers. We showed that in the infrared limit 
the NLIE successfully reproduces the scattering amplitudes involving breathers and 
we gave examples of comparing numerical results from the NLIE to those coming 
from conformal perturbation theory. 

We then proceeded to the case of perturbed superminimal models. Following the 
idea of the introduction of a twist parameter P UHl [16l [TTl [18] , we have shown that 
in the UV limit the twisted NLIE at certain specific values of the twist angle gives 
conformal weights that are consistent with the spectrum of = 1 superminimal 
CFTs. It turned out that similarly to the case of perturbed Virasoro minimal CFTs 
[8] to describe all possible states of the limiting superminimal CFT it is not enough 
to choose a single value of the twist parameter, but a range of values is required 
f l5.6l) . In the IR limit we demonstrated that the twisted NLIE reproduces correctly 
the bootstrap S-matrices of perturbed superminimal models. In the special case of 
SA4{3, 5) + $13 model we checked numerically our NLIEs against the twisted DDV 
equations of ref. [8J and we found perfect agreement. We think that the evidences 
presented in this paper strongly supports the validity of the NLIE description of 
perturbed superminimal models. 

The SSG is just the first in a whole series fractional supersymmetric sine-Gordon 
models [28] that can be seen as perturbations of the conformal models introduced in 
|29j . The full series of coupled NLIEs generalizing the ones presented here has been 
proposed, for the vacuum by C. Dunning [30]. It would be nice to extend Dunning's 
conjecture to excited states as well. To achieve this one has to start from higher 
spin vertex models and along the lines of [TT] and extend the method described in 
[To] for excited states. Having the NLIEs for the excited states two applications can 
be of importance: the quantum group restriction leading to SU(2)-coset theories 
|31j perturbed by $^3^ and the large spin limit where this series should make contact 
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with the N = 2 supersymmetric sine-Gordon model, whose finite size effects are 
of principal importance in the formulation of superstrings propagating in pp-wave 
background [32] • 

Recently in [20] the vacuum state of the SSG with Dirichlet boundary conditions 
has been investigated by means of the NLIE technique. The extension of these NLIEs 
to excited states and to more general boundary conditions would be of interest. 

We hope to return to these issues in the future. 
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Appendix A. 

This appendix is devoted to recall the breather-breather S-matrices of the super- 
symmetric sine-Gordon model [21]. The nth breather has mass 



rrir 



i„ = 2Alsin(^) . (A.l) 
The breather 2-body S-matrix can be written in the form [21] 

where 

^ sinh(g) + t sin(^p7r) sinh(g) + t sin(^p7r) ^ 
^ ^ smh{e) - i smi^pn) sinh(^) - i sin(^p7r) ^ 

'-Zj.i sin2(M^p7r + ) cos2(H±2^j9vr + 2|) 

is the sine-Gordon S-matrix and the supersjmametric factor is 

S^^svi^) = M("'™)(0)G("'™)(^) (A.4) 
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with M("''™')(6') being a 4 x 4 matrix of the form 



/ . sm(Ii£^)+sin(I^) 

^ ' sinh(e) 




^sm(2|E)sin(^) 
cosh(|) 











^ ^ smh(9) ^ 

■ Vsin(I^)sm(ngE) . sin(2g)-sin(H!|Z) 



sinh(|) 





sinh(6>) 







^sin(Iifi)sm(i2f^) 
cosh(|) 




^ ^ . sin(I^)+sin(ir^) 



sinh(6l) 



(A.5) 



and is given by 



(A.6) 



sinh( 



sinh(|) + isin(x7r) 



exp 



q cosh2(f )cosh(7rg) ^ ^ J 

(A.7) 



(A.8) 



Now we restrict our attention to the n = m = 1 case. In this sector there are 4 
eigenvalues of the 2-body breather S-matrix S^^'^^{9). The first pair of eigenvalues 
is of the form: 

a£)(^) = a«(^)g(1'^)(^)4^^)(^), 

where 

The second pair of eigenvalues takes the form: 



(A.9) 



where 



sin^ ffl ^ ^_ 



sinh (^) y cosh^ (f ) 
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Appendix B. 

In this appendix we present the result of the integral (13.21] in case of 5y = 0. Inserting 
(13.61) into (13. 2p and calculating the integral finally one gets: 

/rKc)L^^=,+,^ = if{o - e,,) - Foi9 - e,) - F,{e - e,), (b.i) 

where 

cosh (td£2L) 

Fo(9) = -\n )„J. i +iXGiO). (B.2) 

^ cosh(^) ^ ^ ^ 

and 

^f[(^)-G{e " ' )-G{e )+ Al{Al - At) A'i{Al - A'i) — ^ 

The functions and constants occurring in flB.3p are as follows: 

W{z) = arctanh(e*^"/22) - arctanh(e-^^"/22), (B.4) 

=4D-2 = 2 cosh(^i -62) +4: sin^ j , (B.6) 



A± = VD± /D^, (B.7) 
G{z) = -\nz + ^\n{z*-A^z^ + l). (B.8) 



Appendix C. 



In section Hj in the calculation of the UV conformal weights, the following diloga- 
rithmic sum must be calculated: 



5'o(p) = 2 |l+ {^''2 cos(p)] + L+ [e-3'^2 cos(p)] + 



sin (3p) 
sin (p) 



(C.l) 



where 



X 
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The function Sq{p) has the following simple properties: 



After some algebra one can prove that 



So{p) = — + N,^27r\p-7r 



71 



P 

ttJ 2 



- TT^ I - i NpTT In (4 cos^ p) , 



where 



Np = Np mod 2 A^^ 



[3^1 






- TT. 




-TT- 



(C.3) 

(C.4) 
(C.5) 



and we made the choice of ln(— 1) = ivr. 



Appendix D. 

In this appendix the most important properties of the function Xk{0) will be clari- 
fied. The function Xk{0) is defined as the odd primitive of the kernel 2t: K{0) and 
given by the formula: 

sinh ii- — -\ 

XK{e) = -i\n . XK{e) = -XK{-e) v^ec. (d.i) 

smh (^f + I) 

The branch cuts are chosen to run parallel to the real axis so that Xx(^) be an odd 
real analytic function on the entire complex plane and continuous along the real 
axis. In this case Xk{0) is not periodic anymore with respect to 2m. It is periodic 
only modulo 27r, i.e. the following identity holds: 

Xi^(^ + 27ri)=xxW-27r. 

It follows that the distance between the consecutive cuts is 27ri and the jump of 
Xk{0) is equal to — 27r at each branch cut crossed from below to up. The choice of 
branch cuts is depicted in figure 1. 
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Figure 1: Locations of the branch cuts of Xk{0)- 
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